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This paper presents an approach to place the phasor measurement unit (PMU) 
optimally, which minimizes the setup cost of PMU. This methodology attains 
complete state estimation of the interconnected power networks. An integer linear 
programming (ILP) method is explored for the optimal PMU placement problem. 
It is used to determine the optimal location and minimum number of PMUs 
necessary to make the interconnected power network completely observable. ILP 
may provide many solutions if acquainting buses to zero injection buses are 
unhandled. In the case of more than one solution, a bus observability redundancy 
index and total system observability redundancy index is proposed to find 
the most promising solutions set for redundancy measurement. The proposed 
algorithm is applied to benchmark the optimal PMU placement solutions for 
the IEEE 14-bus, IEEE 30-bus, New England 39-bus, IEEE 118-bus, and NRPG 
246-bus test systems. The obtained results of the proposed approach are 
compared with the existing standard algorithm, and it is observed that the 
proposed approach achieves complete observability of the interconnected power 
network under base-load conditions. 
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1. INTRODUCTION 

Interconnected power network measurement, synchronized on a universal basis, is shifting from 
the laboratory to the utility. PMU-an instrument that uses global positioning satellite (GPS) technology, 
provides novel opportunities for the interconnected power network monitoring, observability, control, 
and protection [1, 2]. In recent years, supervisory control and data acquisition (SC AD A) systems are used 
for surveillance of the network conditions, but it gives unsynchronized measurements dominating 
the inconsistent estimation of the power network states [3]. Furthermore, a scan rate of data (2-4 samples per 
cycle) makes the SCADA system inefficient for measuring the dynamic/transient behavior occurring in 
the power network. These issues can be vanquished by the advent of PMUs, which mitigate this problem by 
using GPS technology and makes precise measurements of the network states [4]. The PMU with a faster 
scan rate (25 samples per second) makes itself acceptable for the observability of the power networks. 
Although the placement of PMU at each bus in the interconnected power network would provide all 
the states of the network system, it is injudicious as PMU and its communication facilities are expensive. 
Thus, an appropriate methodology is necessary for the site selection of PMUs. PMU may measure voltage 
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and current phasors and covers other characteristics such as preservative actions. The goal of the present 
paper is restricted to find the minimum number and the optimal location of PMU for complete observability 
of the interconnected power network states under intact conditions. 

An interconnected power network is said to be fully observable only when all of its states are 
uniquely measured [5-7]. The vigorous research activities on the issues of finding the minimum number 
and optimal locations of PMU have already been published in the open literature. Phadke et al. do pioneering 
work in the field of PMU development and its application in the mid-80s [2, 8]. Some scientists, engineers, 
and researchers believe that the deployment of PMU at each bus will escort to a simplified linear 
state estimator. However, later this problem is resolved in [1], as each PMU can measure not only the bus 
voltage but also the branch current incident to the bus. Hence, proper site selection of PMU can make 
the power network completely observable. In [9 and 10], a novel binary search algorithm is used to find 
the minimum number and optimal locations of PMU for interconnected power network state estimations. 
In [11], the authors used a binary particle swarm optimization technique for finding the optimal locations 
of PMU. In [12], a novel intelligent search-based technique for the placement of PMUs in connected power 
networks while maintaining complete observability is proposed. 

For the optimal placement of PMU, a genetic algorithm-based procedure is developed in [13]. 
In [14, 15], the authors proposed a topologically based three-stage optimal PMU placement (OPP) approach 
for the observability of interconnected power networks. In [16], the authors suggested a novel investment 
decision model for finding the optimal location of PMUs that gives assurance of the complete observability 
of the power grid. An ILP and multi-criteria decision-making based approach is proposed in [17] for placing 
the PMUs in multiple stages over a given period that guarantees fully interconnected power network 
observability even during a line outage or a PMU collapse. In [18, 19], authors used a multi-objective 
biogeography based optimization algorithm for site selection of PMU which makes the power network fully 
observable. In [20, 21], the authors used integer programming to find the minimum number and optimal 
locations of PMU for state estimation. In [22, 23], the problem related to OPP and conventional power flow 
measurements to assure observability during faulted conditions in power networks are considered. In that, 
at the beginning, the methodology is presented as a nonlinear integer programming problem and then 
changed into a similar ILP problem through Boolean suggestions. An integer programming based 
methodology is used by the authors in [24, 25] for the OPP problem in the interconnected power network. 
In [26], the authors explored the consequences of channel capacity of PMUs on their optimal locations to 
ensure that the interconnected power network is fully observable. In [27], the authors proposed a new 
methodology for the OPP problem in a connected power network that is suffering from random component 
outages. In [28], the authors proposed the sum of the variance of the robust estimators to determine 
the PMU placement problem. Also, the placements obtained are further illustrated based on the variance 
of the estimated states. Both the weighted least squares and robust estimators are taken into consideration. 
The OPP problem is evaluated as a binary semidefinite programming model with binary decision variables 
in [29]. Both single PMU and line outage is considered. 

Two deterministic formulations are proposed in [30], which are mixed-integer linear programming 
and nonlinear programming, for solving the OPP problem to achieve complete power network observability. 
The authors have proposed a novel combinatorial formulation for monitoring the complete power grid 
in [31]. A multi-criteria decision support method, analytical hierarchy process, has been used to solve 
the OPP problem. The Pareto approach by nondominated sorting genetic algorithm II is proposed in [32] 
to minimize the PMU placement cost with the current channel selection and the state estimation error. 
In [33, 34], the two-phase branch-and-bound algorithm is proposed for unraveling the OPP problem. 
The main contribution of the presented work is to exclude radial buses (RBs) from the list of potential 
locations for employing a PMU because a PMU located at a radial bus can measure the voltage phasors at 
that bus and only one additional bus which is associated to it. PMU installed at a bus linked with the RB can 
measure the voltage phasor of the radial bus by using the measurement of the current phasor through 
the radial line. Therefore, a PMU is pre-assigned to each bus connected to a radial bus. 

In this paper, we propose an integer linear programming algorithm for finding the minimum number 
and optimal locations of PMUs for the observability of the interconnected power network states. A different 
methodology that is numerical and uses integer programming is conferred. This method enables 
the unchallenging investigation of power network observability for mixed measurement sets. The developed 
criterion has removed any redundancies in PMU placement attained from the proposed algorithm. 
The proposed method is applied on IEEE 14-bus [35], IEEE 30-bus [35], New England (NE) 39-bus [36], 
IEEE 118-bus [35], and northern regional power grid (NRPG) 246-bus [37] test systems. The results obtained 
with the help of the proposed approach have been compared with the existing standard results published in 
the open literature to ensure the efficacy of the proposed method in OPP issues. The remaining paper is 
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divided into the following segments: Section 2 explains the research methodology. In section 3, case studies 
and details of results are given, and in section 4, the conclusion of the paper is presented. 


2. RESEARCH METHOD 

PMU placed at a specific bus in power network precisely measures the voltage phasor of the bus 
and the phasor currents for all branches incident to that bus [9]. Rely on the kind of PMUs used, the number 
of channels used for estimating voltage phasor and current phasor will differ. Therefore, the complete power 
network observability can be achieved by placing a PMU at a specific bus in the system. The main goal 
of the OPP problem is to furnish an observable network by using a minimum number and optimal locations 
of PMUs [38]. For the graphical representation of the OPP, IEEE 14-bus system is shown in figure 1 as an 
example [39]. In the interconnected power network, there is a bus known as a zero-injection bus (ZIB) with 
no generation and no loads connected to it. This ZIB minimizes the number of required PMUs by one 
because it can be merged with one of its adjacent buses. This means that the network topology is modified by 
using the merging process. It needs to redefine the network equations to reflect these changes [40]. 

In IEEE 14 -bus system, three PMUs are located at B2, B 6 , and B9. B7 is only a zero injection bus 
(ZIB). The PMU at B2 not only measure the phase voltage of B2, but also the phase current of Li, L2, L3 
and L4. The phase voltage at Bi, B3, B4 and B5 can be acquired from the line currents by using Ohm’s law. 
After finding phase voltages at Bi, B2, B3, B4 and B5, the phase current of L4, L 6 and L7 are determined. 
Similarly, PMU at B 6 measures phase voltage at B 6 and the phase currents of Lg, L14, L15 and Li6 are also 
measured. Again by using Ohm’s law, the phase voltages at B5, Bn, B12, and B13 can be measured and phase 
current of L19 can be determined. Further, PMU at B9 can measure the phase voltage at B9 and phase currents 
of L11, L12, L13 and L17 and by applying Ohm’s law, the phase voltages at B4, B7, Bio and B14 can be 
evaluated. In addition, the phase current of L9 is also determined. As a result, phase voltages of Bio, Bn, B 13 
and B14 are known, phase currents of Lig and L20 are also estimated. By using phase currents of L9 and L12, 
the phase current of L10 can also be evaluated by using Kirchhoff s current law (KCL). The phase voltage at 
a radial bus, Bg, can now be determined by utilizing the phase voltage at B7 and the phase current of L10. 
Hence, the complete IEEE 14 -bus system is observable by placing PMUs at B2, B6, and B9 and considering 
B7 as a ZIB. 


o 

o 



Figure 1. Optimal placement of PMU for IEEE 14-bus test system [39] 


2.1. ILP based PMU placement method 

The integer programming (IP) is a numerical optimization programming for issues having integer 
variables, and it is the most common method for unraveling the OPP problem. IP is mentioned as integer 
linear programming (ILP) when the constraints and objective function are linear. In an ILP, when some 
variables are integers and other non-integers, then ILP is mentioned as mixed-ILP (MILP). In case 
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the variables are confined within [0, 1] then ILP can be treated as binary- ILP (BILP) technique. Therefore, 
the constraints played a significant role when using the ILP method to unravel the OPP issue. 

2.1.1. OPP problem formulation 

The objective of the OPP is to obtain the minimum number of PMUs needed and their locations for 
achieving full observability of the interconnected power network. Thus, the OPP is formulated as follows: 


n 

mi 

k =1 

s.t. a(Y) > b 


r = Ui y 3 . 


(1) 

( 2 ) 


where, n is the number of buses in interconnected power network for the deployment of PMUs, C/< is the cost of 
PMU set-up at k th bus, Y is the binary decision variable vector having element yu which decides achievability of 
PMU on k th bus and whose entries are defined as in (3). a(Y) are the observability constraints whose entries are 
non-zero if the bus voltage is noticeable w.r.t. the given sets of measurement and zero otherwise. 


J1, if PMU is set - up at k th bus 
[0, otherwise 

The entries in a are as follows: 


a 


k,m 


1, if bus k is linked to bus m 
< 1, if k - m 

0, otherwise 


and 


b = [ 111... Ilf 


(3) 


(4) 


(5) 


The methodology for developing the constraint equation is explained for two attainable conditions 
(1) when there are no conventional measurements or ZIBs and (2) considering ZIBs [39]. The IEEE 14-bus 
test system is taken as an example to describe the above mentioned cases. 

Case 1: A system with no conventional measurements. 

In this case, ZIBs are ignored. For the sake of building the constraint set, the binary connectivity 
matrix a , as defined in (4), is constituted first.st. The a ( nxn ) matrix for the IEEE 14-bus test system is 
given in (6). 


a (nxn ) ~ 


I 1 
1 1 
0 1 
0 1 
1 1 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
-0 0 


0 0 
1 1 
1 1 
1 1 
0 1 
0 0 
0 1 
0 0 
0 1 
0 0 
0 0 
0 0 
0 0 
0 0 


1 0 
1 1 
0 0 
1 0 
1 1 
1 1 
0 0 
0 0 
0 0 
0 0 
0 1 
0 1 
0 1 
0 0 


0 0 
0 0 
0 0 
1 0 
0 0 
0 0 
1 1 
1 1 
1 0 
0 0 
0 0 
0 0 
0 0 
0 0 


0 0 
0 0 
0 0 
1 0 
0 0 
0 0 
1 0 
0 0 
1 1 
1 1 
0 1 
0 0 
0 0 
1 0 


0 0 0 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
111 
0 0 0 
0 0 0 
0 0 0 
10 0 
10 0 
Oil 
Oil 
0 0 1 


Pi 

B 2 

b 3 

b 4 
Bs 
Be 
b 7 
Bs 
b 9 

Bio 
B n 
Pi 2 

Pi 3 
lJ Bl4 


( 6 ) 


The constraints for this condition can be developed as (7). 
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/, 

= y l + y 2 + y5 


>1 

A- 

= >’1 + >’2 + >’2 + 

y 4 y 5 

>1 

/s’ 

= y 2 + y 2 + y 4 


>1 

A- 

= y 2 + y 3 + y 4 + 

y5 + yi + y 9 

, >1 

A 

= y, + y 2 + y 4 + 

y 5 + y 6 

>1 

A : 

= ys + y 6 + y,i + 

■ y t2 + y<3 

>1 

A 

+ 

00 

+ 

+ 

£ 

II 

y 9 

>1 

A 

£ 

+ 

II 


>1 

A 

= y 4 + y 2 + y 9 + 

y 10 + y i 4 

>1 

/io 

= y 9 + yio + yn 


>1 

/.. 

= y 6 + y w + yn 


>1 

/. 2 

= y 6 + y 12 + y i2 


>1 

f\ 3 

= y& + yi 2 + yn 

+ y 14 

>1 

f \14 

= y 9 + >’13 + >’14 


>1 


The operator 4 +’ represents the logical ‘OR’ and the benefit of 1 in the right-hand side of the inequality 
assures that not less than one of the variables present in the sum will be non-zero; e.g. consider the constraint 
related with Bi and B 2 as stated in (8). 


fi = yi + y 2 + y 5 ^ 1 

f 2 = yi + y 2 + y 3 + ^4 + 3^5 ^ 1 

The constraint/i>l signifies that a PMU must be installed either at Bi, B 2 or B 5 to observed Bi. Likewise, 
the constrain/2>1 shows that a PMU must be installed at any one of the Bi, B 2 , B 3 , B 4 or B 5 so that B 2 
becomes observable. 

Case 2: A system with zero injection (ZI) measurements. 

In this case, we contemplate the exceedingly common position where zero injection measurements 
may exist but insufficient to make the interconnected power network completely observable. Again, the IEEE 
14-bus test system, as shown in Figure 2 is considered, where B 7 is a ZIB. It is obvious to identify that if 
the phase voltages at each three out of the four, B 4 , B 7 , Bg and B 9 are acknowledged, then the fourth one can 
be determined by using KCL at B 7 where the total sum of the injected current is known. Therefore, 
the constraints associated with the ZI buses are to adapt with the acquainting buses of theses buses and build 
a set of non-linear constraints and this is concluded as depicted in Figure 2. 



Figure 2. IEEE 14-bus test system with conventional measurements 
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In IEEE 14-bus test system, the constraints associated with ZIB B 7 with its acquainting buses B 4 , Bg, 
and B 9 are altered as given in (9); 

f 4 = 3 2 + 3 3 + 3 4 + 3 5 + 3 7 + 3 9 + /z-Zs'/q ^ 1 

f* = yi + y% + fA-fi-f9 (9) 

fg = y 4 3 7 3 9 3io 3 i 4 f 4 'fi'f% — 1 

The operator V represents the logical ‘AND’ in (9). The explanations for f are clear by utilizing the 
following features of the logical AND (.) and OR (+) operators. A+B=B and A.B-A. Substituting for f from (7) 
in f 4 of (9), yields: 


f 4 = V 2 + 3 3 + 3 4 + 3 5 + 3 7 + 3 9 + fn-fs-fg 

= y 2 + y 3 + y* + y 5 + y 7 + y 9 + (y 4 + y 7 + y 8 + y 9 )-f%-f 9 n m 

= 3 2 + 3 3 + ^4 + V 5 + 3 7 + 3 9 + y 4 -f&-fg + y-j-fz-fg + yz-fs-fg + y 9 -fz'fg ' 1 

= 3 2 + 3 3 + >"4 + V 5 + 3 7 + 3 9 + yz-fz-fg 

In (10), y 4 .f s .f 9 is removed as it is the subset of y 4 . Similarly, y 7 ./ 8 ./ 9 and y 9 .f s .f 9 are also 

eliminated. It is observed that the explanation for f should also indulge an ancillary product term given by 

/ 4 ./ 8 ./ 9 , however this term will be dilapidated. In each simulated condition, this estimation is found 
ineffective in the optimization process. Further, substituting for f in (10), gives: 


f 4 = 3 2 + 3 3 + 3 4 + 3 5 + 3 7 + 3 9 + y%-Uf 9 

= 3 2 + 3 3 + 3 4 + 3 5 + 3 7 + 3 9 + 3^8 * (3^7 + 3 8 )-/ 9 (11) 

= 3 2 + 3 3 + 3 4 + 3 5 + 3 7 + 3 9 + 3 8 -/ 9 

again, substituting for f in (11) yields: 


/ 4 = 3 2 + 3 s + 3 4 + 3 5 + 3 7 + 3 9 + 3 8 -/ 9 

= 3 2 + 3 3 + 3 4 + 3 5 + 3 7 + 3 9 + 3 8 -(3 4 + 37 + 3 9 + 3 10 + 3 14 ) (12) 

= y 2 + y 3 + y 4 + y 5 + yi + y 9 + y 8 .y 10 + 3 8 -3i 4 


similarly, implementing the same logic to other expressions, the constraint set can be written as (13). 

/ 4 = 3 2 + 3 3 + 3 4 + 3 5 + 3 7 + 3 9 + 38-310 + 3 8 -3i 4 ^ 1 

/s = 3 4 + 3 7 + 3 8 + 3 9 (13) 

f 9 = y 4 + 3 7 + 3 9 + 3i 0 + 3 i 4 + 3 2 -3 8 + 3 3 -3 8 + 3 5 -3 8 ^ 1 


It is observed that the constraints analogous to each bus remain unchanged, as shown in (7). But 
the constraint for B7 where the injection is measured is removed from the constraint set as the constraints 
related to ZIBs are incidentally considered by the product terms to add to the constraints related to 
the acquainting buses. The constraint for this condition is given as (14). 


fl 

= Vi + V 2 + v 5 

> i 

fl 

= Vi + V 2 + V 3 + v 4 + Vs 

> i 

fs 

= v 2 + v 3 + v 4 

>i 

f 4 _ 

new = V 2 + V 3 + V 4 + Vs + V 7 + V 9 + Vs-VlO + V 8 -Vl 4 ^ 1 

fs ■ 

= Vi + V 2 + V 4 + Vs + V 6 

> i 

f 6 

= v 5 + v 6 + Vn + yi 2 + y 13 

> i 

fl 

= v 4 + y? + v 8 + v 9 

>i 

fs. 

new = v 4 + y 7 + Vs+y 9 

>i 

fg_. 

= y 4 + v 7 + v 9 + Vio + Vi 4 + y 2 -y 8 + y3-y8+y5 

IV 

fo 

= v 9 + Vio + Vn 

>i 

fii 

= v 6 + Vio + Vn 

> l 

f 12 

= Ve + y 12 + y 13 

>i 

f 13 

= v 6 + y 12 + y 13 + v i4 

> l 

/l 4 

= y 9 + Vis + Vi 4 

> l 


( 14 ) 
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2.1.2. Algorithmic steps to obtain OPP problem 

To find the solution for integer linear/non-linear programming problem, TOMLAB/MINLP with 
MATLAB software package is used [41]. The implementation of the proposed algorithm to obtain OPP for 
an interconnected power network consisting 4 Tl ’ no. of buses using TOMLAB/MINLP is explained as 
follows [42]: 

Step 1. Obtain binary connectivity matrix a(nxn) using (4). 

Step 2. Define b _ l (lower bound on ( y 9 a.y )) and b _ u (upper bound on ( y 9 a.y )). 

Step 3. Explain 'c \ cost of the installed PMUs at k th bus. 

Step 4. Expound’ y _ /, y _u' , bound on variables. 

Step 5. Initialize ’ y ’ vector (y _ 0) , if empty set as 0. 

Step 6. Define' IntVars\ each variable should be integer. This field is explained individually depending on 
the length. Variable indices should be in between [1 n\. 

Step 7. Explain ' flp ', upper bound of f(y ). Only consider solutions < flP - epsilon , default le20. 

Step 8. Interpret the priorities of the integer variables 'VarWeight '. Can be any values, but lower values 
means higher priority. 

Step 9. Assign routine for defining IP problem ' Pr ob — mipAssign{step\ — to — step8 )'. 

Step 10. Display PMU locations 'N PMU ' using tomRun , the universal driver routine. 

2.1.3. Bus observability index (BOI) 

After obtaining the minimum number and optimal locations of PMUs, an expression for BOI is 
given as [39]. We can acknowledge BOI as a performance symbol at all aspects of optimization. 

BOI = a (nxny N PMU (15) 


2.1.4. Total system observability redundancy index (TSORI) 

The TSORI, which is the sum of bus observability for all buses, is given as follows [39]: 


x * PMU 

TSORI = ^ 


i N 

Vnxn) ' iy PMU 


(16) 


where, N PMU is the total optimal number of PMUs obtained from step 10 of Section 3.2, a is the binary 

connectivity matrix which is obtained from (4). As shown in (15) and (16) gives the BOI and TSORI for all 
the possible results of optimal locations of PMUs. 


3. RESULTS AND DISCUSSION 

The proposed approach for finding the minimum number and optimal locations of PMU has been 
applied on standard IEEE 14-bus, IEEE 30-bus, New England 39-bus, IEEE 118-bus and NRPG 246-bus test 
systems. Figure 3, shows the single line diagram of the NRPG 246 bus test system. TOMLAB/MINLP with 
MATLAB software package is used to perform the integer linear programming problem, as described in 
subsection 2.1.2. The technical configuration of the computer is Intel core: I3-2330M (2.2 GHz), L3 Cache: 
3MB, and System Memory: 3GB DDR3 which is used for the simulation purpose. Table 1 shows the number 
and locations of zero injection buses (ZIBs) and radial buses (RBs) for the standard test systems. 
The proposed approach eliminates the RBs from the possible locations of PMUs. As a result, buses next to 
these RBs with higher branch connectivity are preferred as the PMU placement sites [10, 39]. Table 2 shows 
the buses where PMU is essentially placed in order to attain the observability of the RBs. Table 3 shows 
the minimum number of PMUs necessary to make the interconnected power network completely observable 
under usual operating states with computational time. Tables 4 and 5 shows the redundancy measurement for 
the standard test system ignoring and considering ZIBs. 

The bus numbers are designated as B 1 , B 2 , # 3 ,..., B n where n = 1,2,3,..., N. The entries in 
the tables symbolize the number of times a bus phase angle is measured by the optimal locations of PMUs. 
For example, in the case of IEEE 14-bus system, ignoring ZIBs, # 4 is observed three times. B s , B 7 , B 9 are 
observed two times and the rest of the buses are observed once. When considering ZIBs, # 4 and B s are 
observed two times and the remaining buses are observed once. Similarly, for the other test systems, 
redundancy measurement is shown. 


Optimal placement of PMU for complete observability of the interconnected power network... (Rohit Babu) 



142 n 


ISSN: 2252-8792 



Figure. 3 Single line diagram of the NRPG 246 bus test system [37] 


Table 1. Location of zero injection buses (zibs) and radial buses (RBs) 


Test 

Size (buses, 

# of 

ZIBs 

# of 

RBs 

system 

lines) 

ZIBs 


RBs 


14-bus 

14, 20 

01 

b 7 

01 

Bs 

30-bus 

30,41 

06 

Bs^B 9 , B 22, B 25, B 2 7, B 2 q 

03 

Bn, B 13 , B 26 

39-bus 

39, 46 

11 

B±, B 2 , B 3 , B 6 , B 9 , b ±1 , b 13 , B 14 , b 17 , B 19 , B 22 

09 

B30? B 31 , B 32, B 33, B 34, B35, B 3 6 , B 37, 

B33 

Bio, B73, B 87 , Bin, Bii2 ? B lie , Buy 

118-bus 

118,186 

10 

b 5 , b 9 , b 30 , b 37 , b 38 , b 63 , b 64 , b 68 , b 71 , b 81 

07 

246-bus 

246, 376 

58 

B54? b 56 , b 59 , b 61 , b 62 , b 63 , b 69 , b 70 , b 71 , b 72 , 

31 

B 2 , B 4 , b 3 , b 12 , B 30 , b 31 , B 38 , b 41 , B 47 , 




B 731 b 74 , b 75 , b 80 , b 81 , b 86 , b 102 , b 103 , b 104 . 


B 5i» B52? B 53 , B 58 , B 76 , B 77 , B112, B 120 , 




Bio 7 ’ Bi 22 > B 126 , B 129 , B 131 , B 147 , B 154 , B 155 , 


B\ 23 ’ B\ 24 , B i 3 5 , B 149 , B 153 , B 156 , B 159 , 




B\67i b 175 , b 179 , b 180 , b 183 , b 209 , b 2W , b 211 , 

^ 212 ’ B 213 , B 214 , B 215 , B 216 , B 217 , B 2 21, B 2 22, 

^ 226 ’ b 229 , b 230 , b 231 , b 232 , B 2 33, b 234 , b 236 , 
B2371 b 238 , b 239 , b 240 , b 241 , b 243 , B 244 


B\ 72 , B 178 , B i 89 , B 2 q 8 , B 2 24 , B 2 46 
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Table 2. Pre-allocation of PMUs 


Test system Buses where PMU is needed for observability of the radial buses 
_ Ignoring ZIBs _ Considering ZIBs 


14-bus 

b 7 

- 

30-bus 

b 9 , b 12 , b 25 

Bl2 

39-bus 

B'l, B 6 , B 10 , B 19 , B 2 o, B 2 2, B 2 3, B 25 , B 2 9 

B 10? B 20, B 23 , B 25 , B 29 

118-bus 

B g» B 12, B 68 , B 71 , B 86 , B 110 

^12’ B 88 , B 110 

246-bus 

Bs’ B34, B 40 , B 48 , B 54 , B 61 , B 62 , B 63 , B 74 , B 75 , B 80 , 

B 106’ B i22 ? B 125, B 126 , B 132 , B 147 , B 158 , B 168 , B 181 , 

B 185 ’ # 187 » B 191 , B 199 , B 203 , B 219 , B 235 

B 6> B 34 , B 40 , B 48 , B 106 , B 125 , B 132 , B 158 , B 168 , 

B 181 > ^ 185 > ^ 187 ’ B 191, B 199 , B 283 , B 219 , B 235 


Note: ’ entry means none 


Table 3. Optimal locations of PMUs for standard test system with computational time 


Test 

Ignoring ZIBs 

Considering ZIBs 

CPU (s) Elapsed 

system 

# of 

Loc. of PMUs 

# of 

Loc. of PMUs 


time (s) 


PMUs 


PMUs 




14-bus 

04 

B 2 , B 8 , B 7 , B 9 

03 

#2’ #6’ B 9 

0.171 

0.133 

30-bus 

10 

B 2 , B 4 , B 6 , b 9 , b 10 , b 12 , b 15 , B 19 , B 25 , B 2 7 

07 

B 3 , #5, B 10 , B 12 , B 19 , B 23 , B 2 7 

0.140 

0.154 

39-bus 

13 

^2’ B 8 , B 9 , B 10 , B 1± , B 14 , B 17 , B 19 , B 20 , 

B 221 B 23? B 25, B 29 

08 

B 3, B 8 , # 10 , B 16 , B 20 , B 23 , B 25 , # 29 

0.202 

0.193 

118-bus 

32 

B 2 ? ^5’ B 9 , b 12 , b 15 , b 17 , B 2 1, B 23 , B 28 , 

28 

#2’ B± 2 ’ B 15, B 17 , B 21 , B 23 , B 28 , B 34 , 

0.218 

0.238 



^34’ B 37 , b 40 , b 45 , b 49 , b 52 , b 56 , b 62 , b 63 . 


#37’ b 40 , b 45 , b 49 , # 52 , # 62 , b 63 , b 68 , 





^68’ B 74 , b 75 , b 77 , b 80 , b 85 , b 90 , b 94 . 


#7i’ B 75 , # 77 , # 80 , B 85 , # 90 , B 94 , 





Bioi’ B 105, # 110 , ^114 


#ioi> #io5’ #110, B 114 



246-bus 

70 

B 3 , b 6 , b 7 , b 10 , b ±1 , b 18 , b 21 , 

62 

B 3 , b 6 , # 7 , b 10 , b ±1 , b 18 , b 21 , 

0.374 

0.468 



B 22 ’ b 24 , b 29 , b 34 , b 40 , b 44 , b 48 , b 54 , b 56 , 


B 22 ’ B 24, B 29 , B 34 , B 40 , B 44 , B 48 , B 54 , 





B 57 ’ B ei , B 62 , B 65 , B 70 , B 74 , B 75 , B 78 , B 80 , 


#56’ b 57 , b 65 , b 78 , b 80 , b 83 , # 85 , b 91 , 





B 33 , B 85, B 91 , B 92 , B 93 , B 96 , B 98 , B 101 , 


b 92 , b 93 , b 96 , b 98 , b 101 , b 106 , b 109 . 





B 106’ ^109’ ^113’ B 447 , B 121 , B 125 , B 128 , 


Bii3’ B 117, B 12 i, £ 125 , 5 128 , 





B 129’ ^132’ B\3 4 , B 14 q, B 141 , B 142 , B 147 , 


#132’ #134’ #140’ #141’ #142’ 





B 157’ ^158’ ^160’ ^163’ ^168’ ^173> ^181> 


#157? #158’ #160’ #163’ #168’ #173’ 





^185’ B 187, #190’ #191’ #194’ #199’ #201’ 


#181’ #185’ #187’ #190’ B 191 , B 194 , 





#202’ #203’ #215’ B 216 , # 2 19’ B 229 , B 234 , 


b 199 , b 201 , b 202 , b 203 , b 219 , b 229 . 





B 235’ #245 


B 234’ B 235, #245 




Table 4. Redundancy measurement for standard test system with ignoring ZIBs 


Test 

system 

Bus observability index 






























Max. TSORI 

redundancy 

14-bus 

1 , 

1 , 

1 , 

3, 

2, 

1 , 

2, 

1 , 

2, 

1 , 

1 , 

1 , 

1 , 

1 


























3 

19 

30-bus 

1 , 

3, 

1 , 

4, 

1 , 

5, 

1 , 

1 , 

3, 

3, 

1 , 

2, 

1 , 

2, 

2, 

1 , 

1 , 

2, 

1 , 

2, 

1 , 

1 , 

1 , 

1 , 

2, 

1 , 

2, 

1 , 

1 , 

1 










5 

50 

39-bus 

1 , 

2, 

1 , 

1 , 

1 , 

2, 

1 , 

1 , 

1 , 

2, 

3, 

1 , 

2, 

1 , 

1 , 

2, 

1 , 

1 , 

1 , 

2, 

1 , 

2, 

2, 

1 , 

2, 

2, 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 



3 

49 

118-bus 

1 , 

2, 

2, 

1 , 

1 , 

1 , 

1 , 

2, 

1 , 

1 , 

2, 

2, 

1 , 

2, 

2, 

2, 

2, 

1 , 

2, 

1 , 

1 , 

2, 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

2, 

1 , 

1 , 

1 , 

3, 

1 , 

2, 

3 

150 


2, 

1 , 

2, 

1 , 

1 , 

2, 

1 , 

1 , 

1 , 

1 , 

1 , 

2, 

1 , 

1 , 

2, 

1 , 

1 , 

1 , 

1 , 

2, 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

2, 

1 , 

1 , 

2, 

2, 

1 , 

1 , 

1 , 

1 , 

2, 

1 , 

3, 

1 , 




1 , 

2, 

2, 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

2, 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

2, 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

2, 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1, 1 



246-bus 

2, 

1 , 

1 , 

1 , 

1 , 

2, 

2, 

2, 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

3, 

3, 

2, 

1 , 

2, 

2, 

2, 

1 , 

2, 

2, 

1 , 

2, 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

2, 

2, 

1 , 

2, 

1 , 

1 , 

4 

360 


1 , 

1 , 

1 , 

1 , 

4, 

2, 

3, 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

2, 

2, 

3, 

2, 

1 , 

2, 

1 , 

1 , 

2, 

1 , 

1 , 

3, 

1 , 

3, 

1 , 

2, 

1 , 

3, 

1 , 

2, 

1 , 

2, 

1 , 

1 , 

2, 




1 , 

3, 

1 , 

3, 

1 , 

1 , 

1 , 

3, 

1 , 

1 , 

1 , 

2, 

3, 

1 , 

1 , 

1 , 

2, 

2, 

2, 

1 , 

1 , 

1 , 

2, 

1 , 

2, 

2, 

2, 

1 , 

1 , 

1 , 

3, 

2, 

2, 

1 , 

2, 

2, 

1 , 

2, 

1 , 




1 , 

2, 

1 , 

3, 

1 , 

1 , 

1 , 

1 , 

2, 

2, 

1 , 

2, 

2, 

2, 

3, 

1 , 

3, 

1 , 

3, 

1 , 

1 , 

1 , 

2, 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

2, 

1 , 

1 , 

2, 

1 , 

1 , 

1 , 

3, 

1 , 

1 , 




2, 

4, 

1 , 

3, 

1 , 

3, 

1 , 

2, 

1 , 

2, 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

2, 

1 , 

2, 

1 , 

2, 

1 , 

2, 

1 , 

1 , 

1 , 

3, 

1 , 

1 , 

1 , 

1 , 




1 , 

1 , 

2, 

1 , 

2, 

1 , 

2, 

2, 

2, 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

3, 

1 , 

1 , 

1 , 

1 , 

1 , 

2, 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

2, 

2, 

1 , 

2, 

2, 

1 , 




1 , 

1 , 

1 , 

2, 

3, 

1 , 

1 , 

1 , 

1 , 

1 , 

1 , 

1 































Table 5. Redundancy measurement for standard test system with considering ZIBs 


Test 

system 

Bus observability index 































Max. TSORI 

redundancy 

14-bus 

1, 

1, 

1, 

2, 

2, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1 



























2 

16 

30-bus 

1, 

1, 

1, 

2, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1 











2 

31 

39-bus 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

2, 

1, 

2, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1 


2 

41 

118-bus 

1, 

2, 

2, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

2, 

1, 

2, 

2, 

2, 

2, 

1, 

2, 

1, 

1, 

2, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

2, 

1, 

1, 

1, 

3, 

1, 

2, 

2, 

3 

146 


1, 

2, 

1, 

1, 

2, 

1, 

1, 

1, 

1, 

1, 

2, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

2, 

1, 

1, 

2, 

2, 

1, 

1, 

1, 

1, 

2, 

1, 

3, 

1, 

1, 

2, 




2, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

2, 

1, 

1, 

1, 

1, 

1, 

1, 

2, 

1, 

1, 

1, 

1, 

1, 

1, 

2, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1 





246-bus 

2, 

1, 

1, 

1, 

1, 

2, 

2, 

2, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

3, 

3, 

2, 

1, 

2, 

2, 

2, 

1, 

2, 

2, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

2, 

2, 

1, 

2, 

1, 

1, 

1, 

4 

336 


1, 

1, 

1, 

4, 

2, 

3, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

2, 

3, 

2, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

3, 

1, 

3, 

1, 

1, 

1, 

2, 

1, 

1, 

1, 

1, 

1, 

1, 

2, 

1, 

3, 




1, 

3, 

1, 

1, 

1, 

2, 

1, 

1, 

1, 

2, 

2, 

1, 

1, 

1, 

2, 

2, 

2, 

1, 

1, 

1, 

2, 

1, 

2, 

1, 

2, 

1, 

1, 

1, 

3, 

2, 

2, 

1, 

1, 

2, 

1, 

2, 

1, 

1, 

2, 

1, 




3, 

1, 

1, 

1, 

1, 

1 

, 1 

, 1 

, 1 

,2 

,2 

,3 

, 1 

,2 

, 1 

,3 

, 1 

, 1 

, 1 

,2 

, 1 

, 1 

, 1 

, 1 

, 1 

, 1 

, 1 

, 1 

, 1 

, 1 

, 1 

, 1 

, 1 

,2 

, 1 

, 1 

,2 

,4 

, 1 

,3, 




1, 

3, 

1, 

2, 

1, 

2, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

2, 

1, 

2, 

1, 

2, 

1, 

2, 

1, 

1, 

1, 

3, 

1, 

1, 

1, 

1, 

1, 

1, 

2, 

1, 

2, 




1, 

2, 

1 

2, 

1, 

1, 

1 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

2, 

1, 

1, 

1, 

1, 

1, 

2, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

2, 

2, 

1, 

1, 

2, 

1, 

1, 

1, 

1, 

1, 

2, 

1, 




1 , 1 , 1 , 1 , 1,1 
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In the present work, redundancy measurement is also more systematically assigned. For example, in 
the case of IEEE 14-bus the maximum redundancy for a bus is 3 (in case of ignoring ZIBs) and 2 (in case 
of considering ZIBs). The optimal number of PMU obtained using the proposed approach has been compared 
with the existing standard results found in the literature are shown in Table 6. It can be observed that 
the optimal number of PMUs necessary to make the interconnected power network (IEEE 14.30,118, New 
England 39, and NRPG 246-bus system) fully observable are somewhat the same as the results published in 
the open literature. 


Table 6. Comparison of results with the existing standard approach under usual operati ng states 


Test system 

14-bus 

30-bus 

39-bus 

118-bus 

246-bus 

S. Chakrabarti and E. Kyriakides [10] 

4 

10 

13 

- 

- 

T. Venkatesh and T. Jain [12] 

4 

10 

- 

- 

83 

Roy et al. [15] 

4 

10 

13 

32 

- 

S. Mousavian and M. J. Feizollahi [16] 

4 

10 

13 

32 

- 

Sodhi et al. [17] 

7 

- 

26 

- 

135 

B. Xu and A. Abur [21] 

4 

10 

- 

32 

- 

B. Gou [24] 

4 

- 

- 

- 

- 

Proposed approach 

4 

10 

13 

32 

70 


Note: entry means not reported 


4. CONCLUSION 

This paper consummated two goals. Firstly, an ILP based approach for finding optimal locations 
of PMU is explored and secondly, the explored approach is executed and tested on several standard test 
systems. The first step in placing the PMUs is the identification of a candidate location. In a power system, 
there may be certain buses that are strategically important, so that a PMU must be located at each of those 
buses. The remaining buses are made observable by installing a minimum number of additional PMUs. 
The radial buses are excluded from the list of potential locations for placing a PMU. As PMU placed at 
a radial bus can measure the voltage phasors at that bus and only one additional bus that is associated with it. 
Further, PMU placed at the bus connected to the radial bus can measure the voltage phasor of the radial bus 
by using the measurement of the current phasor through the radial line. Therefore, a PMU is pre-assigned to 
each bus connected to a radial bus. Pre-assigning PMUs to certain buses in this manner, the system 
observability is satisfied since no violation constraint exists. Thus, the number of PMUs are the same for 
each power system since the system observability is satisfied with the improvement that no radial buses are 
included in the optimal solutions. The offering of the paper remains in benchmarking the universal OPP 
solutions for a number of extensively used systems and determine the optimal locations of PMUs for 
the system. The obtained results indicate the efficacy of the proposed approach of OPP for interconnected 
power network observability. After comparing with the standard approaches in the earlier published 
literature, the simulation results show that the proposed approach determines the minimum number of PMUs, 
unlike earlier approaches which may find either the same or even higher number of PMUs. 
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